close all;clear;clc

y = @(x)x .* sin(x);
s = @(x)int(f,0,x);
subplot('Position',[0.1,0.1,0.3,0.8])
grid on
fplot(s,[0,4])
xlabel('x')
ylabel('s = \int_0^x(xsinx)dx')